TESTING GRAVITY LAW IN THE SOLAR SYSTEM 
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ABSTRACT. The predictions of General relativity (GR) are in good agreement with observations in 
£f) , the solar system. Nevertheless, unexpected anomalies appeared during the last decades, along with the 

increasing precision of measurements. Those anomalies are present in spacecraft tracking data (Pioneer 
and flyby anomalies) as well as ephemerides. In addition, the whole theory is challenged at galactic and 
cosmic scales with the dark matter and dark energy issues. Finally, the unification in the framework of 
quantum field theories remains an open question, whose solution will certainly lead to modifications of the 
bJT)' theory, even at large distances. As long as those "dark sides" of the universe have no universally accepted 

interpretation nor are they observed through other means than the gravitational anomalies they have been 
designed to cure, these anomalies may as well be interpreted as deviations from GR. In this context, there 
is a strong motivation for improved and more systematic tests of GR inside the solar system, with the 
ON ' aim to bridge the gap between gravity experiments in the solar system and observations at much larger 

scales. We review a family of metric extensions of GR which preserve the equivalence principle but modify 
the coupling between energy and curvature and provide a phcnomcnological framework which generalizes 
the PPN framework and "fifth force" extensions of GR. We briefly discuss some possible observational 
consequences in relation with highly accurate ephemerides. 
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1. TESTS OF GENERAL RELATIVITY (GR) IN THE SOLAR SYSTEM 

The foundations of GR rely on two main pillars. The first one is the equivalence principle which 
states the universality of free fall and gives gravitation its geometric nature. This principle is tested in 
modern experiments at the 10 -12 level, which makes it one of the best tested properties of nature. The 
validity of the equivalence principle has also been tested very accurately in the solar system using Lunar 
Laser Ranging [1] or the Sun-Mars orbit [2], The resulting bound is too small to allow the comparatively 
large anomalies observed in the solar system. Therefore, even if a violation of the equivalence principle is 
nevertheless possible, this strongly indicates that the anomalies, if of gravitational origin, should find an 
explanation in the framework of metric extensions of GR. We focuss on theories which describe gravity 
by a tensor metric field g^. Let us mention at this point that the general static and isotropic metric, 
which can be used as a preliminary description of the solar system, essentially reduces to two functions 
<?oo( r ) and g r r(r) such that 

ds 2 = g^dx^dx" = goo(cdt) 2 + g rr dr 2 (1) 

The Einstein curvature tensor G^ v , built upon the Ricci curvature, has a null covariant divergence 
(Bianchi identity) 

Gfiu = R^-\ g^R , D^G^y = (2) 

This geometrical identity is often put in correlation with the conservation of the stress tensor, D^T^ U = 
0. This remark leads to the second pillar of GR which are the equations relating the curvature of spacetime 
to the energy- momentum content. Those equations allow one to determine the metric tensor from the 
distribution of energy-momentum in spacetime through the Einstein-Hilbert equation which involves a 



unique Newton gravitational constant Gn 
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Note that this form is not imposed by any geometrical argument, so that GR has to be selected out 
from a large variety of metric theories by comparing the predictions drawn from the Einstein-Hilbert 
equation to the results of observations or experiments. When performed in the solar system, the tests 
effectively show a good agreement with the solutions of Q, which means that the metric tensor g M „ has 
a form close to its GR prediction given, for a static point-like mass M by 

GM 

g Q0 = 1 + 20 + 2(j> 2 + ... , -g rr = l-2(f>+ ... with = - (4) 

More precisely, the predictions of the Einstein-Hilbert equation in the solar system are often tested 
in the PPN framework or a Yukawa fifth force framework. The simplest PPN framework is characterized 
by two constant parameters (3 and 7 inserted in the previous Taylor expansion Q 

.900 = 1 + 20 + 2/30 2 + ... , -g„ = 1 - 2 7 + ... (5) 

It turns out that 30 years of tests have put stringent bounds on the parameters /3 and 7 and selected a 
vicinity of GR when analyzed in the PPN framework, namely 7 — 1 = (2.1 ± 2.3) x 10~ 5 and j3 — 1 = 
(—2.5 ± 7.5) x 10 -5 . Concerning 7, the current bound is essentially given by the experiment performed 
through radar ranging of the Cassini probe during its 2002 solar occultation [3] , while the bound on /? is 
obtained via analysis of ephemerides [4] . 

It is important to note that the previous observations have also been tested in other frameworks. For 
example the so-called "fifth force" tests which focus on a possible scale-dependent deviation from the 
gravity force law. Such a deviation, corresponding to an additional massive scalar gravity field, reduces 
to a modification of the Newtonian potential by an additional Yukawa potential 

3oo = [3oo]gr + 20(r)a exp (~t) (6) 

with an amplitude a measured with respect to Newton potential and a range A related to the mass scale 
of the hypothetical new particle which would mediate the "fifth force" . The presence of such a correction 
has been looked for on a large range of distances and it turns out that the Yukawa term is excluded with 
a high accuracy at some ranges, for example a < 10 -10 at A ~ Earth-Moon distance and a < 10 -9 at A ~ 
Sun-Mars distance. These bounds, again deduced from Lunar Laser Ranging and tracking of planetary 
probes, correspond to a remarkable result which approaches the accuracy of equivalence principle tests. 



2. NEW FRAMEWORK OF METRIC EXTENSION OF GR 

We now present a more general framework, which extends the two previously discussed frameworks. 
As already emphasized, the form of the coupling between energy-momentum and curvature can still be 
discussed so that we can generalize Einstein-Hilbert equation ([3]) such that it takes the form of a non local 
response relation between Einstein curvature and the energy-momentum tensor (see [5,6] for details) 

G p „(x) - J dV V "(i - x')T p(T {x') (7) 
As an example, we retrieve GR with the following local expression of the susceptibility x^J" 7 

[ X ,r(x - x')] QR = ^ + 8i5$ 6W(x - x') (8) 

Let us note at this point that generalized response equations naturally arise from radiative corrections 
of GR [5,6]. Radiative corrections may induce modifications of GR not only at high energies, but also at 
large distances [7]. In our simple model of solar system, the stress tensor reads T p<y {x) — S p o5croMc?d^ (x) 
so that the modified Einstein tensor reads 

SG„ u (x) = Mc 2 f cAt 1 S Xflv 00 (x ~ x') , x' = (ct', 0) (9) 



An important feature of this framework is that the Einstein tensor G^, and therefore the Ricci 
curvature do no longer vanish outside the sources, as it is the case in GR. Neverthelesse, because 
GR is in good accordance with observations, it is expected that SG^ should stay small. Within a linear 
approximation, the Bianchi identity can be used to extract the two degrees of freedom 5G^ and SG^ 
which are present : 

»„M = «S5W + «<'>« - ^X^f^ + (10) 

where Tt^uoo an d ""Sjoo are operators of projection on the two sectors of different conformal weights, that 
is to say the traceless and traced part (see [5] for their expression). It has been shown in [6] how these 
two functions SG^ and SG^ are related to a modification Sw of w = ^ In \g r r/goo\ (which is equivalent 
to a modification of the Weyl tensor), and a modification SR of the Ricci scalar R 

5G {0 \r) = -^- i [d 2 r + - r d^5w{r) , SG^(r) = ~^SR(r) (11) 

Generally speaking, in the present framework, the gravitational constant Gat is replaced by a non 
local susceptibility tensor of rank 4 which reduces to only two functions 5G^ and SG^ 1 '. Those two 
functions correspond to two different sectors for gravitation which, in the simplified situation considered 
here, can be exactly matched to the two degree of freedom goo and flVr with the help of (jlip . 

,9oo = [#oo]gr + 5g Q , g rr = [g rr ]GR + Sg rr (12) 

The functions Sgoo(f) and Sg rr (r) can have a general r-dependence. As a pedagological exemple, we 
will focuss in the following on simplified phenomenological models obtained as a Taylor expansion of Sgoo 
and 5g rr in a vicinity of the Solar system {a n r n <C 1 and Xn?" n <C 1) 

<5 5 oo = 2^a„r" , 5g„ = 2 ^ X nT n (13) 

n>0 ti>0 

More elaborated models can be obtained by adding a cutoff to allow the perturbation to start from a 
given distance from the sun. The previous models can be viewed equivalently as distance-dependent PPN 
parameters 

^ (r) " 1 = (^r) £^ n+2 - ^)- 1 = -^E^ ,l+1 ( 14 ) 

^ ' n>0 n>0 



3. PHENOMENOLOGICAL CONSEQUENCES ON A SIMPLIFIED MODEL 

As emphasized in the previous section, modifications of 500 are heavily constrained. Therefore, we 
will consider in the following a modification of the spatial part g rr alone. Then, even if we can deal with 
the general case, we further restrict ourselves to the following simple modification 

5g rr = 2\2r 2 (15) 

This specific model corresponds to a constant Ricci curvature R rr = 8x2, so that \2 represents the 
amount of Ricci curvature outside the source. The geodesic equation is modified, as well as the Shapiro 
time delay, leading to anomalies in spacecraft tracking and ephemerides. For example, a test particle 
(spacecraft or planet) is submitted to the following anomalous coordinate acceleration 

Sa r = 2x2 {rvl - rvj - GM) , 5a s = 4^2 rv r v e (16) 

It is easy to verify that if \2 > 0, a purely radial escape trajectory suffers a positive radial acceleration 
(towards the outside of the solar system) while a circular orbit suffers a negative radial acceleration 
(towards the sun). Moreover, the orthoradial acceleration Sag vanishes for circular orbits and radial 
escape trajectories. The order of magnitude of the anomalous acceleration is given by 



6a - X2GM 



(17) 



This anomalous coordinate acceleration is of course not the observable that is deduced from measure- 
ment. Indeed, one has to take into account the effect on light propagation. The modification ST of the 
one-way light-time T = £2 — t\ from position ri(ii) to position r 2 (i2) can be computed by integrating 
ds 2 = along the path, giving 

cT = R12 + aJ7ihapiro + cST , cST = -^X2i?i2 (r? + r§ + n • r 2 ) (18) 

where fj = \ri\ = [r^GR + tin and -R12 = \r2 — ri\ = [i?i2]cR + SR12 are modified when compared with 
their GR values due to the anomalous coordinate acceleration (TTB|) . The tracking observables can then 
be computed from the knowledge of the two-way light-time (a sum of two one-way light-time connected 
at the level of the spacecraft). For example, up to small relativistic corrections, the Doppler signal 
is obtained by the time derivative of the light-time while the position in the sky (right ascension and 
declination) is obtained through spatial derivatives. 




Figure 1: Geometry of spacecraft tracking. 



An important feature is that the light-time correction naturally produces daily and semi-annual mod- 
ulations, due to the geometry of the triangle sun-station-spacecraft. Those modulations will show up in all 
observables, in particular in ranging. Therefore, correlated anomalies should be observed in cphcmcridcs 
and analyzed within this framework. As a matter of fact, only a careful comparison of observations 
in the outer solar system, within the post-Einsteinian phenomenological framework, could allow one to 
determine whether all gravity tests can be compatible with the anomalies seen in the solar system. This 
comparison will have to account for the presence of the two sectors as well as for their r-dependences. 
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